Is a Classical Language Adequate in Assessing the Detectability 
of the Redshifted 21cm Signal from the Early Universe? 
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The classical radiometer equation is commonly used to calculate the detectability of the 21cm 
emission by diffuse cosmic hydrogen at high redshifts. However, the classical description is only 
valid in the regime where the occupation number of the photons in phase space is much larger than 
unity and they collectively behave as a classical electromagnetic field. At redshifts z < 20, the spin 
temperature of the intergalactic gas is dictated by the radiation from galaxies and the brightness 
temperature of the emitting gas is in the range of mK, independently from the existence of the cosmic 
microwave background. In regions where the observed brightness temperature of the 21cm signal is 
smaller than the observed photon energy, hv — 68(l + 2)~^mK, the occupation number of the signal 
photons is smaller than unity. Nevertheless, the radiometer equation can still be used in this regime 
because the weak signal is accompanied by a flood of foreground photons with a high occupation 
number (involving the synchrotron Galactic emission and the cosmic microwave background). As 
the signal photons are not individually distinguishable, the combined signal+foreground population 
of photons has a high occupation number, thus justifying the use of the radiometer equation. 

PACS numbers: 95.30Dr, 95.55Jz, 98.80-k 
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A black-body radiation field of temperature T and a 
photon occupation number n = (exp{/iz^/fcr} — is 
characterized by fluctuations in the number of photons 
per quantum state with a root-mean-square amplitude 

of H, 



(Arirnis) = \/n{n + 1) = \J v? + ' 



(1) 



At low occupation numbers (n <C 1) or high photon en- 
ergies, the fluctuation amplitude is (Anrms) ~ >A^, as 
expected from Poisson statistics of independent photons. 
This regime is often encountered in optical, UV, or X-ray 
telescopes. However, in the regime where many photons 
occupy the same state In ^ 1), the photons are corre- 
lated and the -n? term dominates over the Poisson term n 
inside the square root of Eq. ([1]). Here, the ensemble of 
photons can be described as a "classical" electromagnetic 
fleld with (Anrms) n- The classical description applies 
to low photon energies hv <ti kT in the Rayleigh-Jeans 
(RJ) regime of the black-body spectrum, and is naturally 
encountered in radio astronomy [l], Q ■ 

The photon number fluctuations induce noise in a de- 
tector which is exposed to the radiation fleld. For a de- 
tector that does not preserve the phase of the incom- 
ing radiation, such as an ideal incoherent bolometer, the 
square of the resulting noise-equivalent-power is given by 
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{An,n,sfhvdiy, (2) 



where A is the collecting area of the detector and is 
the solid angle that it views. For n » 1, Eqs. ([T]) and ([2]) 
can be combined to get the noise level of a measurement 
over a narrow band of frequencies Ah' centered on a fre- 
quency ly = c/A, (AP,.„,s) = 2(^f7/A2)i/2(Ai/)i/2(fcT). 



This result has no reference to Planck's constant h as 
expected for the classical regime, and provides a noise 
power that is simply proportional to the temperature as 
expected for RJ scaling (where for a simple antenna, the 
pre-factor Afi/A^ ~ 1). However if n ^ 1, the Pois- 
son term dominates and yields a very different result 
for the noise power (APrms), which is lower by a fac- 
tor of [hv/kT) exp{— ft,z^/2fcr} than the expression in the 
n ^ 1 regime. In this case, the resulting noise power is 
not simply proportional to the temperature. The preci- 
sion of the measurement is limited by Poisson statistics, 
i.e. by the square-root of the total number of photons 
that are collected by the telescope. Coherent detectors, 
such as interferometers, preserve the phase information 
about the incoming radiation and have a minimum noise 
temperature of hv/k from Heisenberg's uncertainty prin- 
ciple P, 0]. This minimum noise already places them 
on the borderline between the quantum and classical 
regimes. Our subsequent discussion will not be restricted 
to any particular type of detector, but rather focus on the 
unavoidable noise introduced by photon counting statis- 
tics in the incoming radiation field. 

Recently, there has been extensive research on the fea- 
sibility of mapping the three-dimensional distribution of 
diffuse cosmic hydrogen less than a billion years after the 
big bang, through its resonant spin-fiip transition at a 
wavelength of 21cm 0|. Several experiments are cur- 
rently being constructed (such as MWA[2i], L0FAR[23], 
PAPER [sg, 21CMA[3H) and others are being designed 
(SKA[3l|) to detect the theoretically predicted emission 
signal. 

In gauging the detectability of the redshifted 21cm sig- 
nal, the so-called radiometer equation is commonly used 
to assess the sensitivity of an interferometric array of 
low- frequency dipole antennae [l'] . The radiometer equa- 
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tion provides the root-mean-square brightness tempera- 
ture fluctuations (AT^)^/^ of a phased array of dipoles at 
a frequency interval Ai/ centered on a frequency v — c/\, 
in terms of the system temperature Tgyg and the integra- 
tion time iint, 
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where k > 1 is a parameter that describes the overall 
efficiency of the array (and is of order the inverse of the 
aperture filling fraction Q). For an optimized ("filled- 
aperture" ) array of dipoles which are separated from each 
other by ~ A/2, the factor k, is of order unity. The noise 
power provided by the radiometer equation is propor- 
tional to the system temperature as expected from RJ 
scaling, and has no reference to Planck's constant. It is 
therefore valid only in the classical regime. In our dis- 
cussion we will assume that the system temperature is 
dominated by the sky brightness temperature, which in- 
cludes both the signal and its astronomical foregrounds. 
For the sake of generality, we will ignore the technology- 
dependent noise introduced by the detector (be it coher- 
ent or incoherent P, Q), and instead focus on the non- 
reducible effect of the astronomical sky. Our discussion 
applies to an ideal instrument in which the sensitivity is 
dominated by the statistics of photons from the sky and 
not by other (e.g. thermal or instrumental) sources of 
noise within the detector [1, |3| . 

The radiometer equation expresses the simple fact of 
Gaussian statistics that the noise temperature can be 
reduced by the square root of the number of indepen- 
dent samples of this noise [ij. For a frequency band 
Aj/, each random sample lasts a period of time ~ /^v^^ . 
If many photons occupy the same volume of space and 
arrive within the same sample time, then each random 
sampling of the noise would contain many photons. This 
would imply that the noise properties are not dictated 
by individual photons but rather by groups of photons 
arriving within the same time interval ~ Aj/~^ and oc- 
cupying the same region in phase space. Such groups 
of overlapping photon wavepackets can be described as a 
classical electromagnetic wave, whose fluctuation proper- 
ties are not associated with the statistics of uncorrelated 
photons. The photons within each sample are correlated 
with each other because they occupy the same volume in 
phase space, i.e. the same quantum state. 

The radiometer equation is derived using the classi- 
cal language of electromagnetic fields [H. As explained 
above, this language is adequate in the regime where 
many photons occupy each quantum state in phase space. 
This regime, in which photons of a given frequency are 
correlated with each other, is genuinely different from 
the regime often encountered in optical, UV, or X-ray 
observations of astronomical sources, for which Pois- 
son fluctuations of uncorrelated photons determine the 
signal-to-noise ratio 0]. When the occupation number 



is high, there are multiple copies of each photon at a 
given frequency and arrival time, and the detection of 
one photon by a particular dipole does not preclude other 
dipoles from simultaneously detecting another photon at 
the same quantum state. The noise properties are dic- 
tated by the number of random samples, each of duration 
~ lS.v~^ , and not by the number of detected photons. 
This gives no dependence of the noise in Eq. dSI) on the 
telescope diameter (for an optimized array with k of or- 
der unity), in contrast with the dependence encountered 
in optical or X-ray observations, where Poisson statistics 
of uncorrelated photons with different arrival times im- 
plies that the fractional noise amplitude is inversely pro- 
portional to the square root of the number of detected 
photons or the square root of the collecting area of the 
telescope. 

We start our discussion by focusing on the signal pho- 
tons because the foregrounds can in principle be reduced. 
For example. Galactic synchrotron radiation is expected 
to be highly polarized and therefore an appropriate po- 
larization filter can remove its polarized component. (A 
fully polarized foreground could be fully removed while 
losing only 50% of the unpolarized t^,|10,] signal photons.) 
Any such filter should be applied to small patches of the 
sky over which the foreground polarization is nearly uni- 
form. The polarization filter needs to be rotated as a 
function of frequency because of Faraday rotation by the 
interstellar magnetic field, in the intervening ionized gas. 
Such a removal could be of great practical value, since 
a reduction of the foreground brightness by a factor / 
would decrease by a factor of the integration time tint 
required to achieve a particular noise level (AT^)^/'^ in 
Eq. ([3]). For now, let us ignore the foregrounds. 

The number of photons occupying the same quan- 
tum state in phase space around a frequency v de- 
fines the photon occupation number (or the dimension- 
less "phase space density" in bins), n. It is re- 
lated to the spectral intensity, ly (which has units of 
ergs cm^^ s^^ Hz^^ ster^^), through the identity n = 
[c^ /2hv^)Iy^ where h is Planck's constant. The bright- 
ness temperature of a radiation field is defined in radio 
astronomy through the RJ relation, kTf, = (A^/2)/,y (even 
when In is not in the RJ regime), yielding the simple re- 
sult. 
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At the redshifts of interest for the current generation of 
21cm experiments, z < 20, the excitation (spin) temper- 
ature of the neutral intergalactic hydrogen is raised well 
above the temperature of the cosmic microwave back- 
ground (CMB) by the cumulative X-ray and Lya radia- 
tion from galaxies Uj. First, we would like to demon- 
strate that in this regime, the desired signal, i.e. the 
21cm brightness temperature of the gas, is independent 
of the CMB and would have the same value even if the 
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CMB did not exist. This is an important point of prin- 
ciple, since the CMB is characterized by n ^ 1 but the 
signal that we are after could be characterized by n < 1. 
The radiative transfer equation for a spectral line [l^ 
reads. 



dli, 4'{v)hv 
ds 



[712VI21 - ("1-B12 - ^2^21) I„] , (5) 



where ds is a line element, 4'{i') is the line profile function 
normalized by J (f>{v)dv — 1 (with an amplitude of order 
the inverse of the frequency width of the line), subscripts 
1 and 2 denote the lower and upper levels, n denotes the 
number density of atoms at the different levels, and A and 
B are the Einstein coefficients for the transition between 
these levels. We can then make use of the standard rela- 
tions 12l|: B21 = (51/52)^12 and B12 = {g2/gi)A2in/I^, 
where g is the spin degeneracy factor of each state. For 
the 21cm transition, A21 = 2.85x 10~^^s~^ and (72/51 = 3 
[1^ . The relative populations of hydrogen atoms in the 
two spin states defines the so-called spin temperature, Tg, 
through the relation, (71.2/711) = ((72/ffi) exp{— .E/fcTs}, 
where E/k — 68 mK is the transition energy. In the 
regime of interest, E/k is much smaller than the CMB 
temperature Tcmb as well as the spin temperature Tj, 
and so all related exponentials can be expanded to lead- 
ing order. We may also replace Ii, with 2kTb/X^. By 
substituting all the above relations in Eq. we get 

the observed brightness temperature Tjj in terms of the 
optical depth in the 21cm line r ^ 1, 



Th = 



(1 + ^) 
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In an expanding Universe with a uniform hydrogen num- 
ber density A^hi and with a velocity gradient equal to the 
Hubble parameter iJ, 



3 h'^^A2i jVm 
327r £;2 HkT, 
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In the presence of the CMB, (En/k) « Tcmb, and the 
right-hand-side of Eq. ^ is proportional to (T, — Tcmb), 
yielding an emission signal if Ts > Tcmb, as expected. 
But more generally, Eqs. ([6]) and ([7]) imply that as long 
as Ts 3> Tcmb, the magnitude of the emission signal Tf, is 
not dependent on the existence of the CMB (or even the 
particular value of T,, as r cx T-^V Eguations ^ and 
([7]) were already derived before |6l. Il4l. Il5|. 

Allowing for a small-amplitude density fluctuation and 
its corresponding peculiar velocity modifies slightly the 
above expression for the optical depth. For a large-scale 
spherical region of a mean over-density S and a mean 
neutral fraction X^ii at a redshift z, the observed mean 
brightness temperature is 0, Q, , 



Th = 9mK X 



HI 



(8) 



where we have adopted the standard values for the cos- 
mological parameters [l^. The factor of | is obtained 
from angular averaging of (1 -I- cos^ 9), in which the 9- 
dependent term results from the line-of-sight gradient of 
the line-of-sight peculiar velocity in the denominator of r, 
with 9 being the angle between the k- vector of contribut- 
ing Fourier modes and the line-of-sight 15|. Numerical 
simulations of reionization indicate that the brightness 
temperature fluctuations in the 21cm emission of hydro- 
gen at 6 < z < 20 are < 20mK [13, [13]. At redshifts 
z < 6, the observed neutral fraction by mass (mainly re- 
siding in damped Lya absorbers) is Xhi « 3% and so on 
large scales Tf, « 0.3mK(l-H|(5)(l-hz)i/2 ^20^. In compar- 
ison, the observed photon energy of the 21cm transition 
is hv/k = 68(1 + z) imK. 

From Eqs. ^ and ([5]), wc find that n < 1 for cosmo- 
logical regions in which 



Xm < 0.4 
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thus invalidating the use of the classical radiometer equa- 
tion for the signal photons on their own in this regime. 
For example, the detectable 21cm emission after reion- 
ization (smoothed over scales of tens of comoving Mpc) 
20] is entirely in this regime. If the 21cm emission signal 
was the only source of photons on the sky, then its de- 
tectability in the regime where n < 1 should have been 
calculated based on Poisson fluctuations in photon count- 
ing statistics and not the radiometer equation. 

However, the 21cm signal is accompanied by unavoid- 
able foreground photons with n » 1. The foregrounds 
obviously include the CMB, for which 
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hv 



40.1 X (1 + z), (10) 



where Tcmb = 2.73K is the CMB temperature |33| . 
In addition, the Galactic synchrotron foreground has a 
brightness temperature of Tsy„ ~ 140K(7//200MHz)-2.6^ 
yielding a flood of photons with an occupation number 
n « 1.4 X 10'^(7//200MHz)-3-6. Although the polarized 
component of the synchrotron foreground can in princi- 
ple be removed by a polarization filter, the unpolarized 
CMB component cannot be separated from the signal. 
The combination of the signal and residual foregrounds 
results in a population of photons with a high occupation 
number that can be adequately described as a classical 
electromagnetic field. Since the unpolarized signal pho- 
tons are not labeled by a quantum-mechanical tag that 
distinguishes them from the unpolarized CMB photons, 
the detectability of even a weak signal is adequately de- 
scribed by the classical language of the radiometer equa- 
tion. 

We therefore conclude that the use of the radiometer 
equation for gauging the detectability of the 21cm signal 
is justified by the flood of noise photons from the sky that 
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accompanies the seeked-after signal. It is important to 
keep this point in mind when considering; the detectabil- 
ity of other spectral hnes 22, 2^ 24, or continuum 



emission 



Mm fro use cosmic gas, especially at 

wavelengths shorter than a millimeter - where the Wien 
tail of the CMB or other foregrounds are not as bright. 

The Cosmic Background Explorer ( COBE) did already 
operate in the regime where hv/kTcmh > 1, but the high 
temperature of its instruments introduced thermal noise 
3J] that brought it to the regime of n ^ 1 where the 
classical radiometer equation applies. The performance 
of future CMB detectors that will be cooled to tempera- 
tures T < hv/k and will observe the CMB at frequencies 
V > kTcmh/h for which n <C 1, will need to be based 
the quantum-mechanical modification of the classical ra- 
diometer equation where the n term dominates over the 
n'^ term in Eq. ([T]). 
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